Report  No.  69 


Equivalence  of  Voltage  Bias  and 
Geometric  Waveguide  Deisgn  in 
Directional  Couplers 

R.  J.  Feuerstein,  W.  Feng,  J.  C.  Powelson, 
S.  Lin,  L.  Bintz,  and  A.  R.  Mickelson 


January  4,  1996 


Guided  Wave  Optics  Laboratory 
Department  of  Electrical  and  Computer  Engineering 
University  of  Colorado  at  Boulder 
Boulder,  Colorado  80309-0425 


AppsrowS  fcK  p®b®«  \ 


19960126  036 


This  work  was  supported  by  the  Army  Research  Office,  Grant  No,  DAAH04-93-G-0191  and  the  Office 
of  Naval  Research,  Grant  No.  N00014-95-I-0494.  Also,  AT&T  is  providing  a  Doctoral  fellowship  supporting 
J.  Powelson. 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
0MB  No.  0704^0188 


la  REPORT  SECURITY  CLASSIFICATION 

unclassified 


2a.  SECURITY  CLASSIFICATION  AUTHORITY 
_ _  '  DISCASS 


2b.  DECLASSIFICATION /DOWNGRADING  SCHEDULE  „  , , 

N/A 


4  PERFORMING  ORGANIZATION  REPORT  NUMBER(S) 

ECIi:/GW()L/69 


6a.  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL 
University  of  Colorado  (If  appliable) 


6c.  ADDRESS  (City,  StJte,  and  ZIPCode) 

Klectrical  &  Computer  Engineering  Dept. 
Boulder,  CO  80309-0425 


8b.  OFFICE  SYMBOL 
(If  appllcabh) 


lb.  RESTRICTIVE  MARKINGS 

none 


3.  DISTRIBUTION /AVAILABILITY  OF  REPORT 

unrestricted 


5.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 
DOD-ONRN00014-92-J-1 190 


7a.  NAME  OF  MONITORING  ORGANIZATION 
Office  of  Naval  Research 
_ Attn:  Dr.  Arthur  .Iordan,  Code  1114  SI 


7b.  ADDRESS  (City,  State,  and  ZIP  Coda) 

800  N.  Quincy  Avenue 
Arlington,  VA  22217-5000 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 


WORK  UNIT 
ACCESSION  NO. 


8a.  NAME  OF  FUNDING /SPONSORING 
ORGANIZATION 

Office  of  Naval  Research 


8c.  ADDRESS  (C/ty,  State,  and  ZIPCode)  10.  SOURCE  OF  FUNDING  NUMBERS 

«nn  w  PROGRAM  IfPROJECT  iTASK  WORK  U I 

800  N.  Quinry  Avenue  ELEMENT  NO.  NO.  NO  ACCESSIO 

Arlington,  VA  22217-5000 
1 1 .  TITLE  (Include  Security  Clauification) 

(u)  Equivalence  of  Voltage  Bias  and  Geometric  Waveguide  Design  in  Directional  Couplers 


12.  PERSONAL  authOR{S)  Feuerstein,  W.  Feng,  J.  C.  Powelson,  S.  Lin,  L.  Bintz,  and 

_  A.  R.  Mickelson 


13a.  TYPE  OF  REPORT  13b.  TIME  COVERED 

technical  FROM  _  TO 


14.  DATE  OF  REPORT  (Year,  Month,  Day)  15.  PAGE  COUNT 
1996,01,04  12 


*.l 

\ PROJECT 

TASK 

NO 

17. 

COSATI  CODES 

FIELD 

1  GROUP  1 

1  SUB-GROUP 

18.  SUBJECT  TERMS  (Continue  on  reverse  if  necessary  arni  identify  by  block  number) 


19.  ABSTRACT  (Continue  on  reverse  if  necessary  and  identify  by  block  number)  '  . . . 

Using  Poincare  coordinates  we  show  that  any  relative  amplitude  and  phase  between  the 
two  modes  of  a  directional  coupler  can  be  obtained  through  varying  a  single  electrode 
voltage  and  by  selecting  the  correct  length  for  the  coupling  region.  Alternatively, 
curved  waveguide  sections  can  be  employed  to  accomplish  the  same  feat  without  applied 
voltage,  which  we  demonstrate  using  beam  propagation  simulations. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT  21.  ABSTRACT  SECURITY  CLASSIFICATION 

U  UNClASSIFIED/UNLIMITED  □  SAME  AS  RPT.  □  PTIC  USERS _ line- 1  .iss  i  f  lud 

22a.  NAME  OF  RESPONSIBLE  INDIVIDUAL  22b.  TELEPHONE  (Include  Ana  Code)  |  22c.  OFFICE  SYMBOL 

—  R-  Mirkolson  303/492-7539 


DD  Form  1473,  JUN  86  Preeioui  editions  are  obtolete.  SECURITY  CLASSIFICATION  OF  THIS  page 


Equivalence  of  Voltage  Bias  and  Geometric  Waveguide  Design  in 

Directional  Couplers 

R.J.  Feuerstein,  Wei  Feng,  Judith  C.  Powelson,  Sihan  Lin, 

Lou  Bintz,  A.R.  Mickelson 
University  of  Colorado 

Department  of  Electrical  and  Computer  Engineering 
Boulder,  CO  80309-0425 

303-492-7077  robertf@fred.colorado.edu  303-492-3674  (Fax) 

Abstract 

Using  Poincare  coordinates  we  show  that  any  relative  amplitude  and 
phase  between  the  two  modes  of  a  directional  coupler  can  be  obtained 
through  varying  a  single  electrode  voltage  and  by  selecting  the  correct 
length  for  the  coupling  region.  Alternatively,  curved  waveguide  sections 
can  be  employed  to  accomplish  the  same  feat  without  applied  voltage, 
which  we  demonstrate  using  beam  propagation  simulations. 


Acknowledgements.  This  work  was  supported  by  the  Army  Research 
Office,  Grant  No.  DAAH04-93-G-0191  and  the  Office  of  Naval  Research, 
Grant  No.  N00014-95-I-0494.  Also,  AT&T  is  providing  a  Doctoral 
fellowship  supporting  J.  Powelson. 


Feuerstein  et  al.  Equivalence  of  voltage  bias  and  geometric  ... 


Dc  biased,  integrated  optic,  linearized  directional  coupler  optical 
modulators  have  been  extensively  studied^.  They  suffer  from  the 
deficiencies  that  a  large  dc  bias  voltage  with  tight  control  is  required  in 
order  to  maintain  approximately  linear  operation,  and  the  devices  are 
subject  to  dc  drift^  due  to  large  fields  and  low-mobility  ionic  conductivity. 

We  will  describe  a  new  method  to  reach  any  arbitrary  operating 
point  for  integrated  optical  devices  that  use  dc  bias  voltages.  This 
technique  will  greatly  reduce  the  magnitude  or  completely  eliminate  the 
need  for  dc  bias  voltages.  The  technique  may  generally  be  described  as 
the  use  of  waveguide  bends  of  controlled  length  and  curvature  that 
passively  bias  the  modulator  into  another  operating  regime. 

The  propagation  and  coupling  of  optical  modes  in  directional  couplers 
can  be  represented  on  a  Poincare  sphere^.  We  will  use  the  symmetric  and 
antisymmetric  supermodes  of  the  directional  coupler.  The  equations 
governing  the  evolution  of  the  complex  mode  amplitudes  in  a  lossless 
directional  coupler  are^ 

as(z)  =  ejPav  2|aso  cos(bz)  +  j  sin(bz)  -t-  ( i ) 

aa(z)  =  eJ^av  ^|aao  cos(bz)  -  j  sin(bz) 

where  pg  and  pa  are  the  propagation  constants,  Ap=ps-pa,  %=X(V)  is  a 
function  of  the  applied  voltage,  b=((Ap/2)^  and  Pav=(Ps+Pa)/2. 

as(aa)  is  the  complex  amplitude  of  the  symmetric  (antisymmetric)  mode. 

For  a  passive,  lossless  directional  coupler,  as(z)=asoe*^PsZ+(t>0)  and 
aa(z)=aaO®^^^^>  where  ago  and  are  real,  and  is  the  phase  difference  at 
z=0  which  may  be  nonzero  due  to  the  taper  region.  We  assume  adiabatic 
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mode  conversion  in  the  taper  regions  for  the  remainder  of  this  paper;  this 
does  not  qualitatively  change  the  results.  The  Poincare  coordinates^  are 
defined  by 

Si=  Iasl2-laal2  =  Pg-Pg  =  cos(0) 

$2=  2lasaalcos((l))  =  cos((t))sin(e)  (2) 

83=  2lasaalsin((l))  =  sin((l))sin(6) 


where  ^  is  the  phase  difference  between  the  modes,  Ps(a)  is  the  power  in 
the  symmetric  (antisymmetric)  mode  and  {0,4>}  are  standard  spherical 
angular  coordinates  (Fig.  la). 

For  a  two-guide  directional  coupler  we  may  calculate  the  power  in 
each  of  the  two  waveguides  at  the  end  of  the  coupling  region  (z>L)  as 


(3a) 

(3b) 


For  X  =  0,  the  output  state  {81,82,83}  as  a  function  of  L  lies  on  the 
equator  circle  81  =0.  Only  on  this  circle  is  lasl=laal.  Pi  may  vary  from  0  to  1 
while  on  the  equator.  However,  for  some  devices  such  as  linear 
modulators  it  is  necessary  to  operate  off  the  equatorT  In  order  to  vary  81, 
we  need  to  introduce  asymmetry  to  couple  the  symmetric  and 
antisymmetric  modes.  This  can  be  done  using  an  electrooptic  or  other 
material  whose  index  varies  with  an  external  influence.  We  will  now 
analyze  how  to  select  any  point  on  the  sphere  as  the  output  state  of  the 
directional  coupler. 
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The  spherical  notation  defined  in  Eq.  (2)  will  be  used,  where 
<|)=Arg(as/aa),  0<<l)<27t,  and  0=Cos‘^(Ps-Pa),  O<0<  k.  With  asO=aaO=l/i^  and 
(|)0=0,  solving  for  0  using  Eq.  (1)  results  in 

cos(0)  =  1^  -^— I  sin  ^bL)  .  ( 4 ) 

For  X=0,  cos(0)=O  as  expected.  cos(0)  can  cover  the  entire  -1  to  +1  range, 
and  cos(0)=±l  only  when  ApL=m7i/i^  (for  m  odd  integer)  and  XL=±ApL/2. 
The  phase  difference  (}>  at  z=L  is  obtained  from  Eq.  (I)  as 


<^=  tan-l[{(Ap/2+%)/b)}  tan(bL)]  +  tan-l[{(Ap/2-X)/b)}  tan(bL)].  (5) 


For  X=0.  (j)=ApL,  while  in  the  limit  of  large  x,  ({>=0.  Any  point  {0,‘t>}  on  the 
Poincare  sphere  can  be  selected  by  varying  ApL  and  X,  as  we  will  show 
below.  L  and  Ap  can  be  chosen  i  independently  by  the  device  designer. 

Alternatively,  the  evolution  of  the  Poincare  coordinates  as  a  function 
of  L  and  X  is  given  by, 
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where  Sio=Si(z=0).  This  shows  that  points  on  the  sphere  are  reached 
through  successive  rotations  about  the  S3,  Si,  and  S3  axes.  With 
{Sio, 820,830}  ={0,1,0},  Eq.  (6)  reduces  to 
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2  ^ 

{Si,  S2,  S3}  =  {  (1  -  cos(2bL)),  (I)  +  (^)  cos(2bL),  (^)  sin(2bL)  } 

!iZl3 


(7) 


We  find  from  Eq.'s  (3a)  and  (7)  that 
Pi  =l-{^)  sin^bL). 


(8) 


From  Eq.  8  we  see  that  dPi/dx  =  0  when  x=0,  which  is  why  a  bias 
point  with  is  necessary  for  a  linear  modulator.  Let  us  now  examine  the 
power  Pi  for  ApL=7i.  When  %=0,  Pi=0  as  expected.  Pi  reaches  its  maximum 
of  one  when  bL=m7t,  m  integer.  Thus  by  varying  X,  Pi  can  be  set  to  any 
desired  value  between  0  and  1.  For  2X/AP=0.8,  Pi=0.5  and  the  slope  of  Pi 
versus  X  is  approximately  linear,  a  good  operating  point  for  a  linear 
modulator.  Then  {Si,  S2,  S3}  =  {0.8,0,-0.6}  which  is  indicated  by  the  *  in  Fig. 
la. 

Another  representation  for  evolution  of  the  output  state  {81,82,83}  is 
shown  in  Fig.  2.  The  starting  points  are  all  on  the  rim  of  the  circle  where 
X=0.  As  X  grows  large  the  trajectory  is  an  ever  tightening  spiral  about  an 
axis  which  approaches  the  S3  axis. 

With  {810,820,830}  ={0,1,0},  Eq.  (7)  may  be  solved  for  the  X  and  L 
needed  to  reach  the  output  state  {81,82,83}. 
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This  shows  explicitly  how  any  point  on  the  Poincare  sphere  can  be  reached 
by  a  suitable  choice  of  %  and  L. 

As  shown  above,  setting  the  operating  point  of  a  symmetric 
directional  coupler  device  off  the  equator  requires  a  nonzero  %.  We  will 
now  describe  a  method  to  reach  any  point  on  the  Poincare  sphere  with 
little  or  no  dc  bias,  by  using  a  waveguide  configuration  which  breaks  the 
symmetry  along  x,  thus  providing  coupling  between  the  symmetric  and 
antisymmetric  modes. 

The  technique  may  generally  be  described  as  the  use  of  waveguide 
bends  of  controlled  length  and  curvature  to  passively  bias  the  modulator. 
Figure  3  inset  shows  the  waveguide  design  with  sine  bends  for  low  loss^  as 
defined  by 

=  (11) 

where  {x,z}=  {0,0}  is  at  the  inflection  point  of  the  individual  waveguide 
bends  and  Lb  and  h  are  defined  in  the  Fig.  3  inset. 

The  beam  propagation  method^-^  was  used  to  simulate  mode 
propagation,  with  light  starting  in  waveguide  1.  Figure  3  shows  how 
Si=cos(0)  varies  with  the  bias  section  length  Lb-  The  graphs  are  flipped 
about  the  Si=0  axis  when  the  direction  of  the  bend  is  up  instead  of  down. 

It  may  be  that  the  limiting  values  for  Si  of  ±1  are  not  achievable  using  this 
technique.  For  example,  in  the  case  of  h=50  |im  and  Lb=2760  [im,  we  find 
Si =-0.996,  (or  +0.996  for  the  up  bend  case).  For  practical  devices  this 
range  should  be  more  than  adequate. 
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Once  the  bias  section  ends,  the  phase  difference  <|)t  between  the  two 
modes  may  not  be  the  desired  value  (|)bias-  A  symmetric  section  of  parallel 
waveguides  with  length 

L(|>=  (‘I’bias'^l’t)  /(Ps"Pa)>  <l’bias>‘l’t  (12) 

(*l>bias+  27C-<t>t)  /(Ps"Pa)»  <l’bias<‘l>t 

changes  (|)  to  the  necessary  value.  This  demonstrates  that  almost  any  bias 
point  on  the  Poincare  sphere  can  be  reached  with  only  passive  waveguide 
design.  Any  displacement  in  the  actual  bias  point  from  the  design 

point  due  to  fabrication  variations  can  be  corrected  with  a  small  dc  bias 
voltage  in  the  L<|)  section. 

We  have  shown  how  a  dc  bias  voltage  and  careful  selection  of  the 
device  length  on  an  electrooptic  directional  coupler  can  set  the  operating 
point  of  the  device  anywhere  on  the  Poincare  sphere.  We  have  also  shown 
that  this  can  be  accomplished  through  the  use  of  geometrical  waveguide 
design  with  zero,  or  a  greatly  reduced  bias  voltage.  This  relieves  the 
requirement  for  large  dc  bias  voltages  and  their  attendant  problems  in 
electrooptic  devices. 
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Fi£ure  Cantinns 

Figure  1.  (a)  Definition  of  Poincare  coordinates,  (b)  Geometry  of  coupler. 

Figure  2.  Trajectories  on  the  Poincar6  sphere  as  a  function  of  XL  (from 
0-10)  projected  onto  the  S2-S3  plane.  Various  ApL  values  are  indicated. 
The  circle  represents  possible  states  when  X  =  0.  The  lower  axis 
represents  the  power  Pi  in  waveguide  1. 


Figure  3.  Si  vs.  transition  length  for  the  passively  biased  waveguides 
shown  in  the  inset.  The  parameters  used  are:  n  =  1.5,  An  =  0.005,  G  =  5.1 
pm,  TE  mode,  X  =  1.3  pm,  h  =  50  pm  and  h  =  75  pm. 
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Figure  1.  (a)  Definition  of  Poincare  sphere  parameters,  (b)  Geometry 
of  coupler. 
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Figure  2.  Projection  of  trajectories  on  Poincare  sphere  to  the  S2-S3  plane  as 
%  increases  from  0-10.  Shown  for  various  ApL  values  as  indicated.  The  circle 
is  %  =  0.  Also  shown  is  an  axis  for  P^  . 
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Figure  3.  Plot  of  versus  transition  length.  The  parameters  used 
are:  n  =  1.5,  An  =  0.005,  G  =  5.1  pm,  TE  mode,  X  =  1.3  ^im,  h  =  50 
|im  and  h  =  75  p-m.  Inset  shows  geometry. 
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